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ABSTBACT 


Measurements of thermal Ions are Influenced by factors such as 
spacecraft potential, velocity, angle of attack, and sheath size. A 
theory Is presented for the response of an Instrument which accepts Ions 
only within a small angle of Incidence from normal. Although a more 
general theory Is available and forms the basis of this one, the small- 
angle restriction allows a simpler formulation which does not depend on 
sheath size. Furthermore, practical Instruments are easily designed 
around this restriction. Laboratory tests verify that such Instruments 
respond as expected and they Illustrate how design details Influence per- 
turbations from the Ideal response characteristics. 
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INTRODUCTION 


In this report Is proposed a design for a spacecraft instrument which 
detects charged particles. Its function is like that of a Faraday cup yet 
its design alleviates many of the problems associated with Faraday cups. 

The instrument described here provides a measure of ion current entering 
an orifice where only those ions having energies above some threshold are 
counted. In th^s way it differs from a class of Instruments , curved-plate 
electrostatic analyzers, which respond to particles within a band of 
energies. Though the theory pertains to either electrons or ions, the 
implementation described here is for ions. 

The traditional Faraday cup has generally consisted of an orifice 
and a collector plate separated by several biased grids. Current reaching 
the collector is detected by an electrometer as a programmed sequence of 
voltages is applied to the grids. Problems in the use of these Instruments 
have been numerous and are detailed in Appendix I. Of particular interest 
for ion collection are sensitivity to photons, grid-wire spacing effects, 
Mach effects, acceptance cone, sheaths, and response time. 

Generally Mach effects and sheaths will be present and should be 
accounted for in the theory of the device; a way of doing this has been 
presented by VRilpple, Warnock, and Winkler (Ref. 1). Included in their 
calculations is the effect of a wide-angle acceptance cone i^lch is 
typical for a Faraday cup. 

This report considers a special case of the aforementioned work, that 
where ions are received only through a narrow acceptance angle as measured 
from the normal to the face of the Instrument. The restriction in angle 
has several beneficial effects, not only simplifying the theoretical 
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formulation but also allowing a simple Instrument design to be Implemented. 

If Ions can approach a biased grid from various angles, the criteria for 
passing through the grid depends on angle as well as Ion energy. Restric- 
tion to near-normal Incidence thus removes one complicating factor In the 
theory. Another feature of the general theory, the dependence of Instrument 
response on sheath size. Is also made unnecessary by the small-angle 
restriction. These simplifications allow replacing the numerical quadratures 
of the general theory with closed-form solutions. The proposed design 
allows for eliminating grids, for separating photons and Ions to render 
the Instrument Insensitive to photons, and for achieving high sensitivity 
with a particle counter which Is used In place of an electrometer. 

In this report the general formulation of the problem Is modified 
for the small-angle case, the Integrals are evaluated, possible applications 
are discussed, and laboratory tests of prototype Instruments are presented. 
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THEORETICAL FORMULATION 


The equations and notation of Whipple, Warnock, and Winkler (Ref. 1) 
are taken as the starting point for this development. Their constraints 
on the range of Integration are Included and then further restrictions are 
Imposed by the specifications on angle of Incidence. With these restric- 
tions, upper and lower bounds on the Integrals can be established. When 
the limits are sufficiently close to each other, a simplified Integral 
can be used for practical computations. 

Summary of Work by Whipple. Wamock, and Winkler 
The current to an Ion measuring Instrument Is 

1 - AJ J (1) 

o 

where A Is the area of the orifice, Is the random current density given 
by 

J - ne(kT/27rm)^^^ (2) 

o 

and where 

J ■ exp(-Vg-M^) II F dx dz 

Allowed 

X,Z 

• exp(-Z-X+BE^^^ cos o^) sin a„) (3) 

The function F Is assumed to be unity In the referenced work and the 

notations, B-2M cosy and D*2Mslny, are newly Introduced here. Other 

symbols are as previously defined and are summarized In the symbol table. 

Symbols are defined with reference to the system shown In Fig. 1. It Is 

especially Important to note the condition of energy conservation 

E - X + Z + V„ (4) 

s 

where all quantities are energies normalized by the factor kT, X and Z are 
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Fig. 1. A spacecraft model which shows a coordinate 
reference system used for defining symbols. 
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tangential and normal kinetic energies of a particle at the orifice of 

the instrument, and V is the normalized energy derived from spacecraft 

8 

potential. Mach number is M and sheath size is imbedded in a as well as 
the integration limits. 

A spherical conducti.re spacecraft of radius R is assumed to have a 

sheath defined in terms of the parameter L > y/R where y is Debye length 

and L varies from 0 to <». This parameter relates the angler of incidence 

measured at the orifice a and measured outside the sheath a . Note 

0 

2 

that tan ■> X/Z. Specifically it may be shown that, for L ■ 0, 




sin a - 
00 



sin a 


E sin ■ X 
00 


(5) 


E cos a » Z+V 
00 3 


and, for L ■ 



a 


00 O 


sin a^, large E 
, small E 


( 6 ) 


VIhen Vg < 0 (attractive to ions) and ~ sin a^, the ratio of cx^/a^ 
is as shown in Fig. 2 . Thus for some I'Oint X,Z one may infer that 


sin a < sin 

00 — 00 L*0 


(7) 


cos > cos 


The allowed range of integration is subject to constraints. First 

one must have E ^ 0 or X+Z -V^. This boundary is Identified by the 

parameter in Fig. 3. In addition the requirement that trajectories 

not start elsewhere on the spacecraft places a further constraint shown 

by the various L parameters. When is positive one simcly requires that 

X > 0 and Z > 0, and in any case the repeller bias condition is that Z > V , 

8 

where V is the normalized repeller voltage. 

O 
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Fig. 2. The ratio of a /a for various L values and small a's. 
° 00 o 



-Vs X 


Fig. 3. Allowed ranges of integration for Eq. 3 for various 
L values. 
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Function Bounds 


When L>0 the Modified Bessel function simplifies as follows: 

I„(D-^ sin ajjj^ . I^(dA) 

Consequently for any L, by use of Eq. 7, 

1 < I (D»^ sin a ) < I (Di/Sc) 

— o flo'_0 

It is convenient to use a series representation 

_4„2 _6„3 

In similar fashion one may write 
B cos a 


®|L=0 

Also by definitions 

B < b/F +WZ 
— s 

where W = 1+U and U H (X/Z) 


= B/V +Z . 
8 


( 8 ) 


(9) 


( 10 ) 


( 11 ) 


( 12 ) 


max 


Thus for any L, one may write 

bA +Z < B E^^^ cos a < Bi'V -WZ 


(13) 


s — 

Small-angle Restriction 

When the Instrument is designed so that ions can enter only at some 

small angle a ^ a^, then the range of integration is modified as shown in 

Fig. 4. This is the same as Fig. 3 with th'. additional restriction that 

X ^ UZ. A conceivable design to limit the angle is to place two circular 

orifices in a row such that only particles passing through both would be 

detected. Fig. 5 shows this arrangement where the smaller of the two 

2 

orifices fixes area A = irD2/4 and the larger of the two fixes U: 




(14) 


This result is applicable if case the parameter F in 

Eq. 3 is unity. 
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Fig. 4. Allowed ranges of integration for Eq. 3 when a constraint 
is placed on the angle of incidence. 
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Fig. 5. Trajectory of particle passing through upper and 
lower circular orifices having diameters and 
^2 respectively. 
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A more practical system uses two orifices of equal diameter as shown 

2 

In Fig. 6. Here A » 'irD /4 as before but U differs: 

U = (|)^ (15) 

Most Importantly the function F Is unity only for a«0 and otherwise It Is 
less than unity. 

For nonnormal Incidence, the effective area Is less than that of a 

single hole. Figure 7 shows the geometry on which Is based a calculation 

of the reduction factor F. 

The parameters 3 and C are found from 

C = S(tana )/2 (16) 

o 

and 

B = arccos (2C/D). (17) 

Then the area A* can be expressed as 

A' = B D^/4 - tanB (18) 



Fig. 6. Trajectory of particle passing through equal 
circular orifices having diameter D. 
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Fig. 7. Computation of the effective area A" which Is 

seen by particles having an angle of Incidence a^. 


10 


Ihe foreshortened area A" Is given by 
A" = 2A' cos 

which is the effective area seen by particles at the angle 

F is found by comparing A" with A cos and it is found to be 

_ 6D^/2 - 2C ^tanB 

F - 2 

ttD /4 


or 


F - -|(2B-8in(26)). 


Equations 15-17 combine to show that 
B“arccosT=Tr/2 - arcslnl 


where T = 


Series expansions in T yield 


B 


»/2. + + 

^^^23245 


euid 


•) 


sin 2B = 2T»^5^ = 2T(1 - -|t^ - -|t^ - ...)■ 


Consequently the final expression for F is series form is 


( 19 ) 

The function 


( 20 ) 


( 21 ) 


( 22 ) 

(23) 


F - 1 - 4) T + (|f) + (^) T= + (^) l' 

1 9 

+ { — - — ^ T +*»* 0<T<1 

^57.6ir^ ^ ^ 1 ^ 1 

4 2 

where the coefficients ~ » 'JJf * **' conveniently designated Fj^, 

F^, *’*. Though this form seems unwleldly compared with Eq. 20 it is 
well suited when the Integral of Eq. 3 is evaluated. The function F is 
Illustrated in Fig. 8. 

Basic Integral 

The inequalities of Eqs. 9, 13 may be used to set upper and lower 
bounds on the Integral of Eq. 3, which is rewritten here in a form 
adaptable to identifying the bounds: 


(24) 


J 


exp(-Vg-M^) J dZ exp(-Z+Bv^ZW+Vg) 
UZ 

I dX f(^| exp(-X) I^(D*^) 


( 25 ) 


o 
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F 



Fig. 8. The function F(T) which represents the decrease 

In aperture area associated with nonnormal Incidence. 


The lower limit V Is to be the largest of -V , 0, and V where V Is the 

normalized potential of the spacecraft relative to the space plasma and 

V Is the normalized grid potential relative to the spacecraft. Note that 
8 

In practice a grid Is not used but rather a biased drift tube. 

For special case where -V > V , - V >0, and L 0, the range 

8 g 8 

specified for the Integral, Eq. 25, will be Incorrect as Fig. 9 Illustrates. 
The contribution AJ from the double-hatched region will be bounded as 
follows : 

0 

•exp(V^ + Ib|/^) I (d/^) (W )^/2 

9 8 0 8 8 

This correction, which depends on the second power of U, Is usually found to 
be negligible. Consequently one may choose the limits from V to <» as If 
L“0. 


£ AJ < exp(-V^-M^) I I dX dZ. 


(27) 
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'Vs X 


Fig. 9. The region excluded from Che approximate Integral 
formulation given by Eq. 25. 


Aside from Che possible contribution of AJ one may then specify bounds 
on Eq. 25 as shown In Table 1. The second case shown Is more simple to 
evaluate than the others and It can be used to calculate the Instrument 

response when the bounds are acceptably close together. The explanation for 

the existence of the bounds comes from Che Inequalities previously noted. 
Though W has been defined as l-HJ so as to set the upper bound In Eq. 13, 
one may obtain the lower bound from W«l. In like manner the upper bound In 

Eq. 9 corresponds to use of D as it Is defined. Setting D*0 produces the 

lower bound. 

In this section the Integral of Eq. 3 has been transformed Into Eq. 25 
which Is a form appropriate for analzying response of an Instrument having 
a restricted angle of Incidence. Table 1 summarizes several special cases 
of Interest and Eq. 27 provides a correction to the upper bound. The function 
F shown in Eq. 24 Is appropriate for the case of two equal aperatures, and 
It Is used In the next section which describes the evaluation of the Integral. 
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TABLE 1. Choices of Parameters which Yield Upper and Lower 
Bounds . 


B 

D 

W 

VALUE OF 
INTEGRAL 

> 0 

Itself 

1 + U 

/ 

max 

0 

1 

J . 
min 

< 0 

Itself 

1 

J* 

max 

0 

1 + u 

“^mln 


^Increase by AJ from Eq. 27 to obtain upper 
bound. 
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EVALUATION OF THE INTEGRAL 


The Integral Eq. 25 contains two variables of Integration, X and Z. 
Integration over X Is done first and the result Is expressed as a power 
series In U, which Is the maximum value of X/Z. When U Is small, as Is the 
case for the prototype systems described In a later section, then the first 
term of the series often suffices. The Integral over Z is done analytically 
without further approximations being made. 

X- Integral Ion 

A portion of Eq. 25 may be Identified as the X- integral: 


J 


X 



exp(-X) I^(D.^) 


(28) 


0 ^ ' 

If the series expansions (see Eq. 10, 24) for I^, F, and exp(-X) are multiplied 
together and (28) Is Integrated term by term, then the results may be collected 
by powers of U: 


>2,d2 


(29) 


+ (UZ)^[f - 1] + 1 Fi + y F 3 + f ^5 •••] 

+ (UZ)^ [— — + — 1 f— + — F + — F + — F •••] 

Note that the F^'s are the coefficients of Eq. 24 for the case of two equal 
apertures and the F^'s are all zero If one aperture is large compared with 
the other. The integral may be written simply as 


Jx (UZ)^Aj (30) 

where the A^'s are given In Table 2. Values up to F^ have been retained In 
evaluating the A^'s because the F-series is slowly convergent. 

Z- Integration 

When from Eq. 30 Is substituted into Eq. 25, the expression for J 

becomes ^ 

00 

J - exp(-Vg-M^) j dZ exp(-Z+Bi^r+i«)^|^ (UZ)^Aj^. (31) 

V 
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TABLE 2. Values of 


One 

aperture 

much 

larger 

than 

other 

^2 

1 

'(i‘. 1) (i) 

'64 4 2^ '3^ 



0.249 

Equal 

A- 

2 

(7 - 1) (0.061) 

apertures 

2 




4 -2 1 
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A change of variable la Instituted through the definition: 

2 V 

J + z - -H« 

4 W 8 

This may be solved for Z explicitly: 

2 V 

Z - K + jj- + B/wh2 

From this the differential dZ Is found to be 
dZ - (2H -f B«^) dH. 


(32) 


(33) 

(34) 


In these equations B may be either negative or positive depending on whether 

the Mach angle y exceeds 90° or not. In terms of the new variable H the 

expression for J Is ^ 

V 2 f 

J - exp(-V + + ^) J dH exp(-H^)(2H+Bi^) 


(H^ + »^_-j& + hB^)^ 


where the limit G Is 


2 V 

G - (-^ if v-bA^+mv) 


1/2 


(35) 


(36) 


This limit of Integration may be negative for certain conditions as shr>vn 
In Fig. 10. Always the limit V will be subject to the conditions 


V > 0 , V > - 

The expression for J contalrs Integrals of the form. 




exp(-H ) dH 


(37) 


(38) 


several of which are shown In Table 3. The Integrals are used as needed 

In developing the expressions for J which follow. 

Equation 35 for J may be expressed In terms of J^'s as follows: 

V 7 

J - exp(-V^ + -ji - + ^) A^U^ (39) 

Typically only !■! will be of interest yet expressions for both axul J 2 
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G 

1 



Fig. 10. Evaluaclon of the lower limit G Illustrated for the 

caae W*l, V —1, B«+l. 

8 — 

are given below. The general form for J. la 

( 2 V 

- J dH exp(-H^)(2H+Bv^)(H^+HB*^ * ^ ~ 

G 

or in simpler notation 
00 

- |2dH exp(-H^)(H+a)(H^+2aH+b)^ 

G 

where a ■ B»^/2 and b ■ a^ -V /W. 

8 

Conaequently one finds that 

- (2G^ + 6aG + 2 + 4a^ + 2b) 1^^ 

+ (3a -f 2ab) I 

0 

J 2 • (2G^ + 10 ac3 +(4 + 16a^ + 4b )G^ 

+ (30a/2 + 8a^ + 12ab) G 
+ (4 + 16a^ + 4b + 2b^ + 8a^b)] 

V fl5a/2 + 4a^ + 2ab^ + 6ab] I 


Ti^LE 3. Integrals 


^0 

t/if er£c(G)/2 

^1 

exp(-G^)/2 

^2 

GI, + 1/2 
i 0 

^3 

(G^ + 1) 1^ 

I4 

(G^ + 3G/2) + 31^/4 

'5 

(G^ + 2G^ + 2) 


g”‘^ I^ + (n-1) 


The integration is thus complete in a general form i^icb appears some- 
what cumbersome but which is adapted in rather simple form to numerous cases 
described next. The most general solution is formed by combining Eq. 42, 

43 for and J 2 with the expression for J given as Eq. 39. The coef- 
flents come from Table 2, the limit G from Eq. 36, and the integrals Ij^, 

1 from Table 3. 
o 
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SPECIAL CASES AND APPLICATIONS 


The generiLl formulation presented in the previous section may be 
interpreted in numerous ways. Several special cases are described here: 
spherical probe, unbiased stationary sphere, biased spheres with 0° and 
90° Mach angles, and the general equations divested of the parameters 
used tc estimate bounds. Then typical parameters are summarized for various 
measurements which have been reported in the literature. 

Sperical Probe 

The equations for a spherical probe may be derived for the special case 
that V - 0, l.e., there is no dlscrimli atlon in terms of the energies of 

O 

particles entering the orifice of the Instrument. Also it is assumed that 

the probe is stationary in the plasma, M=B=D=0. Furthermore F is assumed 

to be unity. If the probe is positively biased (V > 0) then V = 0, and 

s 

for this repulsive mode of operation: 

= Vg/W, Aj^=l, A2 = - I (44) 

Also: 

Jl = = 4Ij (45) 

Consequently, with reference to Eqs. 1, 39^ 

I = AJ exp(-V ) (U - •••) (46) 

o s 

This is the usual expression for current collected by a repulsive probe 

2 

except for the factor (U-U • • • ) which accounts for the restricted angle of 
incidence. 

For the case where the probe attracts (V < 0) , one uses V = -V and 

s s 

2 

G = V /W - V . The Integral then becomes 
S 8 X 

= exp(-G^)(l-V^) (47) 

and if J2 is ignored then I becomes 

I = AJ^ (l-Vg)U. (48) 
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Thl% aside from the factor U, is the equation for an attracting spherical 
probe . 

Unbiased Stationary Sphere 


Another relatively simple case has V = 0 and V >0 where again, 

s g 

M=B»DsO and The response for V < 0 Is the same as for V =0. For 

8 8 

this case: 


Aj^=l, A 2 = - 2 

Jf = (1+V ) exp(-V ) (49) 

J„ = (V ^+2V +2)exp(-v ^ 

2 8 8 o 

The current Is found to be 

I = AJ exp(-V )[U(1+V ) - U^(V ^/2 + V +!)•••] (50) 

o 8 8 8 8 

Biased Sphere with Zero Mach Angle 

For this case it is assumed that ^ -1 making the sphere attractive, 
that M= 1/2, and that 7*0. Consequently B=l and D=0. It is also assumed 
for illustration that U = 0.01 amd that F=l. The lower bound for the 
integral corresponds to W=1 such that 


Aj^=l, A 2 = - 1/2, a= 1/2, b= 5/4 

g 2 = - 3/4 + V - /-1+V 
8 8 


(51) 


Using these values one may show that 

= exp(-G^) (G^ + 3G/2 + 11/4) 
+ lli^ erfc(G)/8. 


(52) 


A similar calculation may be done for W=1.01 which corresponds to the upper 
bound. Results for both values of W are shown in Fig. 11 which includes 
only the term. Inclusion of the J 2 term at the point = 1 yields (with 
G= - 1/2) 

I = AJ (5.458 U - 9.047 U^ •••) 
o 
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2 


i 




Fig. 11. Bounds on J/U for U=0.01, B=l, D=0, and F=l. Only 
is included in the calculation. 



2 

from which it is seen that the J 2 term containing U represents a correction 

of less than 2 percent. The correction AJ from Eq. 27 applies only if 

Vg < 1 and it has an upper bound of 

AJ/U < 4.3x10"^ (53) 

\dilch is about 0.1 percent. If one is willing to Ignore the various uncer^ 

tainties amounting to at most 3%, then he can use the case for W=1 without 

J 2 to write an expression for curi-'ni as a function of grid bias: 

I = AJ U J, (54) 

o 1 

Here and G are as defined by Eqs. 51 and 52. 

The numbers and results shown in this case are considered to be a typical 
application though of course the equations cease to be useful if, for some 
combination of parameters, the bounds are not sufficiently tight. The 
general equations as formulated in this report allow one to estimate the 
possible error in the use of a simpler equation such as (54) . 
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Biased Sphere with 90° Mach Angle 


Parameters are the same here as for the previous case except that 

B*0 and D=l. Consequently W«-l, G*A -1, A.» 1, a*0, and b*l. The para- 

8 

meter A 2 depends on D which will be 1 for the upper bound and 0 for the 
lower bound. Thus A 2 may be -3/8 or -1/2, respectively. Consequently the 
Instrument response Is 

I = AJ^ exp(-G^-l/4) [U(G^+2) -U^(G^/2 + 2G^ + 5/2^***] (55) 

If V = - V which Is the condition for maximum I, then G=0 and 
8 s 

I = AJ^ exp(-l/4) [2U - 5U^/2 •••] (56) 

The above expressions yield the lower bound with D=0, the upper bound Is 

2 2 
found by reducing the U term by the factor 3/4. lAien the U term Is 

Ignored, the result becomes 

I = AJ^ exp (3/4 - Vg) U (V^ + 1) (57) 

Note that for the special case B^O, the value of W Is Immaterial. The same 
result Is obtained for either W=1 or W=l+U. 

Simplified Equations 

When one has Justified that the bounds are sufficiently tight by 
considering the various cases of Table 1, then It Is expedient to work with 
the response In simplified form. To do this one assumes W-1 and J2=J2***=0. 
Consequently one has 

G^ = B^/4 + V + V - bA +V (58) 

s s 

and the constant A^ may have either of two values, 1 or 0.249, depending 
on the choices of aperture diameters. Then 

I = AJ^ exp(-M^+B^/4) Aj^ U (59) 

where 

J » exp(-G^)(G^+3BG/2 + 1 -V +3B^/4) (60) 

X 8 

+ v¥ erfc(G) (3B/4 + B^/8 - BV /2) 

s 
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In these expressions the value of B may range from -i2M to -2M as the 

Mach angle changes. As previously given, V is the largest of 0, -V , and 

8 

2 

V . The value of G will be negative if V < B /4 - V and B > 0. 
g ® s 

Possible Uses 

Numerous measurements have been made where the proposed instrument and 
its theoretical interpretation would have been applicable. Some of these are 
summarized here. 

Whipple, Warnock, and Winkler 11 were able to adjust the para- 

meters of ion temperature, spacecraft potential, and sheath size until 
their theoretical predictions matched the data of OGO 3. Ion density was 
also found from the normalization factor J^. The equations developed here 
can be used similarly except that sheath size is not needed nor can it 
be Inferred from the data. 

Whipple, Wamock, and Winkler [Ref. 1] refer to OGO 3 measurements 
where ion temperature varied from 4000 to 20,000 K and the corresponding 
factor kT/e from 0.34 to 1.7 ev. For satellite potentials to -5V, the 
normalized potential -V could be as large as 16 in the plasma sphere but 
less than 1 beyond the plasmapause. Mach number M was typically less than 
unity and substantially less beyond the plasmapause for hydrogen. The 
spread on the bounds of the Integral Increases with B and with ~Vg> Typical 
percentage uncertainties in J are given in Table 4. Obviously one would 
try to keep spacecraft potential from being much larger than the ion 
temperature (in eV) if small error bars are sought. 

Harris [Ret. 2 ] advocates the measurement of cold ion density in the 
plasma trough as an Important diagnostic tool and, in a related paper, 
Chappell [Ref. 3] describes detached regions of the plasma sphere which 
have densities several hundred times the usual trough density of about 


24 


6 3 

10 ion/m . An Ion sensor as described in this report could provide 
Information about the cold (le, thermal) Ions In the plasma sphere and 
beyond . 

Sensitivity analysis can be based upon the simplified form given by 

—6 2 

Eqs. 59, 60. Assume that A ■> 0.785x10 m^, A^ » 0.249, and U - 0.0062, 
where these values correspond to a prototype which has been laboratory 

“16 

tested. Random current density at 4000°K Is, for hydrogen, « 3.68x10 n(A) 

_3 

or, with a particle counter, 2300n (Hz) tdiere n Is In m . For an M of 

0.5, a B of 1, and V = - V = 10, 

’ ’ g s ’ 

I = 0.0380 A J = 69xl0“^n (Hz), 
o 

Because the background counting rate Is typically 0.2 Hz, then a density 
4 -3 

of 1.4x10 m will be detectable at a rate of about 1 Hz. At the other 

extreme, counting rates as high as 10^ Hz are possible without pulse over- 

9 -3 

lap, this corresponding to a density greater than 10 m . The sensitivity 
of an Instrument can be set through the choosing of the aperture area and 
the separation between apertures; a pair of Instruments can cover comple- 
mentary ranges. 
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TABLE 4 


Typical Percentage Uncertainties In J Values 


-V 

f, 

V 

.8 

M 

B 

u 

Percent Difference 

J -J . 
max min 

Percent 
Error if J 
Neglected 

10 

10 

1 

2 

0.01 

10. 

*3.1 

10 

10 

1 

2 . 

0.0062 

6.2 

1.9 

10 

10 

0.5 

1 

0.0062 

3.7 

1.8 

10 

10 

0.5 

0 

0.0062 

0.0 

1.5 

1 

1 

0.5 

1 

0.0062 

0.7 

0.5 

1 

4 

0.5 

1 

0.0062 

0.8 

0.8 


DESIGN AND TESTING OF A PROTOTYPE 


The preceding theoretical development provides a basis for the design 
of an Instrument to measure thermal Ion density and temperature. The key 
requirements are that the Instrument accept Ions from a small viewing cone, 
that It provide an Internal retarding potential for rejecting Ions having 
energies below an adjustable threshold, and that the entrance orifice be 
biased more negatively than the local plasma potential unless some other 
means exists for measuring spacecraft potential. Then from the data and 
with a calibrating measurement of sensitivity, one may calculate the Ion 
density, the Ion temperature, and the spacecraft potential. 

The retarding potential In the Instrximent should be stepped through 
a sequence of values ranging from zero to, say, 5 V more than the sum of any 
negative bias on the Instrument Itself plus the maximum negative excursion 
In spacecraft potential which Is anticipated. However, for resolution one 
wants perhaps several steps per volt. Thus a compromise arises Involving 
resolution, range of spacecraft potential, and telemetry rate. The use- 
fulness of the Instrument can be extended to times when negative charging 
would otherwise disrupt operation by operating an electron emitter In such 
a way as to stabilize spacecraft potential. 

Instrument Design 

A prototype Instrument has been constructed as shown In Fig. 12 to 
demonstrate the feasibility of detecting Ions subject to the criteria of 
small viewing cone and energy selection. The Instrument was enclosed In 
a grounded stainless steel can, which was supported on pivots. The can 
could be rotated about an ends passing through the entrance orifice and 
lying In the plane of the orifice. 
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Fig. 12. Diagram of a prototype ion sensor showing an Ion 

trajectory, orifices, repeller, and an Ion detector 
which viay be biased up to 3kV. 
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Two holes of 1-mm dlam and separated by 12.7 mm defined the orifice 
area A and maximum angle of incidence or alternatively the parameter U, 
which was 0.0062. After Ions passed through the orifices, they entered 
a drift tube maintained at potential v if their energy was sufficient. 
Otherwise they were collected on the back side of the orifice assembly. 

Ions passing down the drift tube were drawn out through a port to 
the negatively biased electron multiplier, each ion having a certain prob- 
ability of producing an avalanche which could be counted. Note that low 
energy ions will have nearly equal trajectories and energies at the funnel 
so that they will have nearly equal detection efficiencies. Tlie exception 
is for ions of high atomic mass which are of little Interest for the 
intended application. Because of the geometry, the funnel will not intercept 
electrons streaming through the orifices nor will it attract them, though 
it may collect ions generated by electrons entering the Instrument. Also 
the funnel is placed so as not to Intercept photons entering the Instrument. 
Vents in the case were shielded to prevent particles from entering other 
than through the Intended orifices. Particles and photons passing through 
the drift tube were collected on the end plate though with some reflections 
occurring. A more carefully designed dump should ultimately be used. 

The original Instrument was fabricated with the two entrance holes 
erroneously misaligned by about one hole diameter. Measurements were 
taken with that configuration and then the outer hole was bored out and 
covered with an allgnable diaphragm having the same hole diameter as before. 
Additional measurements were made for the new configuration (modification of 
Fig. 12) which had a hole separation S>15.9 mm and U*0.0040. 

Test Sources 

Two test sources could be mounted Interchangeably with their apertures 
30 cm above the face of the ion sensor. One was a commercial Lyman«a 
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source producing an estimated 1.5x10^^ photons cm ^s ^ at the face of the 
instrument. The other was a heated tungsten filament which could be biased 
negatively to produce electrons and positively to produce Ions by contact 
Ionization. The charged particle source, shown In Fig. 13, had a grounded 
cover and a focusing cup which was connected to the more negative side of 
the filament. The spread of particle energies from the source depended 
upon the length of the hot tungsten wire and the associated potential drop 
along Its length. 

The ion source operated with ambient pressures from 10 ^ to 4x10 ^ 
torr when air was being leaked into the system. No attempt was made to 
measure what constituents were being Ionized, but presumably predominant were 
O 2 and H 2 O which had the lowest Ionization potentials of the common 
constituents of air. The source would not function with helium which has 
a much higher ionization potential. The ion source was operated with bias 



Fig. 13. Crojs section of ion source. 
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vrltages from 6 to 30 volts. Stable and consistent operation was obtained 

-13 

et the higher voltages and typical current density at the sensor was 10 
A/ cm^ . 

Response to Ions 

Measurements were made with Ion beams having different energies, 
different energy spreads, and different angular orientations. 

The response at normal Incidence as a function of repeller voltage 
was measured for two Ion beams having different accelerating potentials. 
Results are shown In Fig. 14 where the low counting rates are attributed 
to an assumed loss of source Intensity when accelerating voltage Is low. 

Figure 15 shows a steeper transition as well as the stronger signal associated 
with a source bias of 30 V. The Ion energy spread depends however, not on 
the source bias, but on the length of the filament, which was reduced to 
6 mm for the test shown In Fig. 15. This cutoff Is characteristic of a 
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Fig. 14. Average counts for time periods shown when the Ion source was at 
5.9V and lOV. 
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DRIFT TUBE POTENTIAL (Vg ) 

Fig. 15. Ion sensor response with an Ion source bias 
voltage of 30V. The slope of a semllog plot 
corresponds to an Ion t«nperature of O.lSeV. 
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Boltzmann factor, exp(-eV/kT), where T is 0.15oV. Consequently distributions 
such as diose in Fig. 14 represent spreads in particle energy, not instru- 
ment response characteristics. 

When the instrument is rotated, the response should vary according to 
the function shown in Fig. 8. Keeping only the linear portion one has 

F - 1 - 4T/7T (61) 

The graph in Fig. 16 compares the measured response with the straight line 
approximation for the original system with the misaligned holes. An 
asymmetrical response is noted. However, this distortion was eliminated 
by aligning the holes. 

Another type of problem with the response is Illustrated in Fig. 17 
which is for the system with the holes aligned. Here the angle corresponding 
to peak response la found to be a function of repeller voltage, which 
Ideally should not be the case. The explanatici; for this effect is pre- 
sumably related to the asymmetric placement of the collection funnel. Though 
this asymmetry may be a problem for some types of measurement, it should 
not affect the response to isotropic sources such as the thermal plasma of 
interest. 

Close examination of che data shown in Fig. 17 for a repeller voltage 
of 30.7 indicates still another type of distortion. The curve is redrawn 
in Fig. 18 to a different scale. The response on either side of tho max- 
imum is proportionally less relative to the peak than for the other cases. 

Both of the previous problems are probably related to the nature of 
the potential contours in the drift tube. Ions must pass through a saddle 
point and ions with barely enough energy to pass may be reflected if they 
are not directed toward the optimum point. Thus ions entering at an 
angle will have a lesser probability of passing consistent with the data for 
30.7 volts. Ions having energies considerably in excess of the minimum 
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ROTATION (DEGREES) 

Fig* 17. Instrument response when the ion source is bi/«sed at 30V. 
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the ion source Is biased at 30V. 


will pass regardless of the exact nature of the saddle point as Is the 
case for the data at 0 volts. 

Response to Electrons and Photons 

The response to electrons depends upon entering electrons producing 
Ions by collision with background gas or the walls. The response Is 
strongly dependent on cleanliness and It decreases with time. Also careful 
allgxuaent reduces the response. No quantitatively valid data was recorded 
because of the variability of the response but counting rates eventually 
dropped below 1 Hz with the Instrument facing the source. 

The response to photons was high, with direct llliunlnatlon giving 
rise to a count rate of 800 Hz. This Is equivalent to approximately 3000 Hz 
In sunlight. However, the angle of response was limited to about +5**. 

With care In designing the dump, one can reduce this rate considerably 
though how much Is not known. The Instrument was sensitive to leakage 
through the vents; counting rates of 5 Hz were obtained with the orifice 
covered. A small side lobe In the response occurred 14° off center when the 
entrance holes were misaligned. 

Alternate Design 

A different drift tube design has been described by Ross (Ref. 4) and 
It Is Illustrated In Fig. 19. After Ions pass through the double apertures, 
they pass through a saddle-point potential which Is controlled by the bias 
on the short drift tube between ground planes. The Ions then exit from 
the chamber through a third aperture. Ross showed that If the third 
aperture Is too small, the electrostatic focussing produced by the drift 
tube causes anomolous response. The dashed line In Fig. 20 Is that 
response which can be avoided by making the aperture large enough. For 
this system the energy of particles which can pass through the saddle point 
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Fig. 19. Schematic of an alternate drift tube design. 

Is some fraction of the potential applied to the tube, that fraction being 
specified In Table 5. 


TABLE 5 


Ratio of Drift-Tube 
Radius to Length 

Ratio of Saddle-Point 
Potential to Applied Potential 


-2 

00 

ir 

2.0 

0.146 

1.0 

0.232 

0.5 

0.721 

0.25 

0.974 

0.20 

0.992 


38 



150 



ONODaS/SiNflOO 


39 


Fig. 20. Normal and anoi 
19.3eV and the 


CONCLUSION 


An instrument for measuring temperature and density of thermal Ions 
above the plasma pause has been designed, built, and lab-tested, as 
detailed In the preceding chapters of this report. The design Is an 
application of the general theory for Ion measurements In the presence 
of sheaths and Mach effects [Ref. 1]. The predominating feature of 
the design Is that particles having large angles of Incidence are excluded; 
in practice the angle is limited typically to 4®. As a consequence of 
this feature, sheath size has little significance and theoretical approx- 
imations yield relatively simple expressions for Instrument response. When 
prototypes were tested with monoenergetlc sources, various abnormalities 
were Identified and associated with design flaws for which corrective 
procedures are prescribed. 

Results of the theoretical development are summarized by Eqs. 1, 36, 

39, 42, and 43 which yield I, the current being measured, as a function 

of the many different parameters which may Influence It. Table 2 provides 

the coefficients A^ used in these equations. Seldom Is It necessary to 

use this full set of equations. Eqs. 58, 59, and 60 together represent 

a simplified set of equations which yield I In closed form and In terms of 

the parameters B, M, V , V, U, J and A. Typical error bars are given 

s o 

In Table 4. 

Because the angle of Incidence Is limited, the energy discrimination 
process Is easy to Implement. A simple biased drift tube provides resolution 
as low as 0.15eV as Illustrated by the data In Fig. 15. In fact the 
energy spread associated with potential drop along the source filament Is 
easily noted In Fig. 14. 
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The theoretically predicted angular response function In the shape 
of the triangle was observed In numerous tests. Various minor anonolles 
In that shape can be attributed to design features which may be optimized. 

A saddle-point discriminator originally described by Ross [Ref. 4] Isolates 
that region from the collection region such that the designer would have 
more control over the anomolles. 



APPENDIX I 


SOURCES OF ERROR IN THE 
USE OF FARADAY CUPS 

These devices, known as traps, Faraday cups or retarding potential 
analyzers, have been used extensively both above and below the plasma pause, but 
the results have often been found to be In error. In retrospect. Cures are 
known for most of the problems, but the construction and Interpretation 
of a fully debugged device Is a major task. 

The device generally consists of a collector plate separated from 
the plasma by numerous biased grids. It Is flush mounted on the surface 
of the spacecraft. Current to the collector Is monitored with an elec- 
trometer as a sequence of programmed voltages Is applied to the grids. 

Problems arising In the use of these types of Instruments are Itemized 
below. However, some of these are applicable to other kinds of particle 
detectors as well. 

1. Internal photoemlsslon. When sunlight strikes the collector. It 
produces photoelectrons which can be returned to the collector 
by a suppressor grid. However, reflected or oblique light 
strikes the back side of the suppressor and photoelectrons emitted 
there are attracted to the collector. Thus a spurious response 
occurs. Serbu and Maler (Ref. 5) Ignored the data which was taken 
when their trap faced the sun (within 70°). Bridge et al. (Ref. 6) 
modulated their particle beam with a driven grid so that the 
desired current could be separated from the steady photocurrent 
with a tuned amplifier. 

2. Photoemlsslon and secondary emission from the spacecraft. These 
electrons are extremely difficult to distinguish from natural 
electrons, especially with a cup. They are not very Important 
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at low altitudes where natural densities exceed, say, 10 cm , 
but they are dominant above the plasma pause. Rosenbauer (Ref. 7) 
considers a method of separating the two types of electrons with 
a curved plate analyzer. 

3. Secondary emission from suppressor grid and collector. The 
suppressor Inhibits loss of secondaries from the collector, but 
It produces Its own secondaries when struck by the Incoming beam. 
Most of these will migrate outward and In any case, their signif- 
icance can be kept low If the suppressor Is highly transparent. 

4. Secondary emission from walls of cup. This can be minimized by 
careful design. 

5. Ion Impact on retarding grid. When electrons are being collected. 
Ions may Impact on the retarding grid so as to produce secondaries 
which are attracted to the collector. This problem occurs only 
below the plasma pause, and It may be eliminated by proper design. 
Knudsen and Harris (Ref. 8) and Blnsack (Ref. 9) have described 
the effects. Blnsack used It as a plasma pause detector. 

6. Grid wire spacing. Temperature calculations should be corrected 
to account for the error In assuming that a grid Is equivalent to 
an equlpotentlal plane. Hansen (Ref. 10) describes the correction. 

7. B Fields. These are generally Ignored and with Justification 
though Wrenn and Helkklla (Ref. 11} have studied the effects on 
photoelectrons at low altitudes. 

8. Ion species. If dlffer?-nt t;i)es of Ions are present (0+, He-t-, 

H4- at low altitude) then special computations are necessary 
(Ref. 12). See also Item 9 which is related to this effect. 
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9. 


Mach effects and wakes. These are Important especially at low 
ali^ltude where Ion temperature Is low, average Ion mass Is high, 
and spacecraft velocity Is high. Above the plasma pause, hydrogen 
ion velocities at 1 keV are comfortably in excess of spacecraft 
velocity but problems still occur for low energy ions. Whipple 
et al. (Ref. 1) include this effect in their theory. Sagalyn and 
Smiddy (Ref. 13) treat the problem as do Serbu and Maier (Ref. 5). 

10. Aperture and acceptance cone. These depend on grid transparency, 
grid spacing, particle type (possibly), and particle energy. 

They can be measured and/or calculated. See Bonettl et. al. (Ref. 14). 

11. Backscatter. Electrons approaching the collector may backscatter 
with sufficient energy to escape. The fraction is typically from 
10 to 30%, and it is larger for high Z metals (Ref. 15). 

12. Potential and sheath effects. The response of a trap depends on 
spacecraft potential and sheath size. These effects may be 
treated by iterating to a solution with potential and sheath 
thickness as parameters. The calculations developed by Whipple 
et al. (Ref. 1) involve numerical quadratures. The effects dis- 
appear when potential goes to zero, when tenq>erature becomes 
large, or when sheath size goes to zero. 

13. When thermal densities are low, the signal associated with thermal 
ions may be obscured by the constant current associated with the 
high energy flux. 

14. The sensor must be mounted and spacecraft attitude controlled so 
that the bearing of the trap does not change during the settling 
time . Otherwise data taken in directional plasma flows may be 
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averaged in an unpredictable manner. The ability to reorient 
the trap la desirable. 

15. Conducting surfaces should have uniform work functions If low 

energy measurements are to be accurate. Gold plating and carbon 
have been used to achieve uniformity. 
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